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Abstract. Boltzniann-Gibbs measures generated by logarithmically correlated 
random potentials arc multifractal. We investigate the abrupt change (" pre-freezing" ) 
of multifractality exponents extracted from the averaged moments of the measure 
- the so-called inverse participation ratios. The pre-freezing can be identified with 
termination of the disorder-averaged multifractality spectrum. Naive replica limit 
employed to study a onc-dimcnsional variant of the model is shown to break down at 
the pre-freezing point. Further insights are possible when employing zero-dimensional 
and infinite-dimensional versions of the problem. In particular, the latter version 
allows one to identify the pattern of the replica symmetry breaking responsible for the 
pre-freezing phenomenon. 
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1. Introduction 

Investigations of multifractal measures of diverse origin is a very active field of research in 
various branches of physics for several decades [1]. To set the notations, consider a certain 
(e.g. hypercubic) lattice of linear extent L in A^— dimensional space, with M ~ standing 
for the total number of sites in the lattice. The measures of interest are usually defined 
via weights < pj < 1 associated with every lattice site i = 1,2, ... ,M and normalized as 
Si^i Pi — 1- One can imagine a few different spatial arrangements of weights pi accross the 
lattice sites. In the case of simply extended measures the weights are of similar magnitude 
at each lattice site, the normalisation condition then implying the scaling pi ~ in the 

large-M limit. As a generalisation of the above example one can imagine the non-zero weights 
Pi supported evenly on a fractal subset of lattice sites of effective dimension < N^j < N. 
In the limiting case of N^j = we then deal with localized measures characterized by the 
weights Pi essentially different from zero only inside one or few blobs of finite total volume. In 
such a situation weights stay finite even when M — > oo, that is pi = 0{M°). Finally, in the 
most interesting case of multifractal measures the weights scale differently at different sites: 
Pi ~ M~"» [§|. The full set of exponents < Oj < oo can be conveniently characterized by the 
density p(a) = '}2i=i ^i'^ ~ '^i) whose scaling behaviour in the large-M limit is expected to 
be nontrivial: p{a) ~ M^'^'^\ with the (convex) function /(a) known in this context as the 
multifractality spectrum. Note that the total number m(a) = p{a) da of sites of the lattice 
characterised by the scaling exponents < a must satisfy the inequality m(a) > M-l'^"^ > 1, 
hence f{a) > 0. The condition f{a) = defines generically the minimal q;_ and maximal a+ 
threshold values of the exponents which can be observed in a typical realisation of disorder. 
Note that the constraint Pi < 1 implies a_ > 0. 

An alternative, frequently more practical way of describing multifractality is via the set of 
exponents Tq characterising the large-M behaviour of the so-called inverse participation ratios 
(IPR's) Pq which are simply the moments of the corresponding measure: 
M „ 

Pi = T.Pi= M-''y{a)da. (1) 
1=1 

Substituting in the above definition the relation p(a) ~ M-^*^°) one can evaluate the integral in 
the large-M limit by saddle-point method. One then finds the relation between Tq and /(a) 
given by the Legendre transform: 

Pg~M~^9, Tq = qa-f{a) where Q = ^ ■ (2) 

The relation is valid as long as the q— integral is dominated by the saddle point. It is easy to 
see however that at large enough |g| the integral should be dominated rather by the vicinity 
of the thresholds a± resulting in linear behaviour of the exponents with q, i.e. Tq = qa±. 

The above description is valid for multifractal measures of any nature. In recent years 
important insights were obtained for disorder-generated multifractality, see [2j and also [3J for 
a comprehensive discussion in the context of Anderson localisation transitions, and [1] for an 
example related to polymers in disordered media. One of specific features of multifractality 
in the presence of disorder is a possibility of existence of two different sets of exponents, 
Tq versus fg, governing the scaling behaviour of typical Pq versus disorder averaged IPR's, 

§ Usually one defines exponents via the relation pi ^ j^-^ui j g_ ^^iq reference to linear scale L 
instead of the total number of sites M ^ . We however find it more convenient to get rid of trivial 
spatial dimension factor TV, and concentrate only on essential parameter behaviour. 
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< >~ M^'^t. Here and henceforth the brackets stand for the averaging over different 
reahsations of the disorder. The difference is related to a possibihty of disorder-averaged 
moments to be dominated by exponentially rare configurations. A related aspect of the 
problem is that the " annealed" multifractality spectrum recovered from the disorder-averaged 
multifractal exponents fq via the Legendre transform ([1]) can be negative: /(a) < 0. Indeed, 
those values reflect events which are exponentially rare [5] and need exponentially many 
realisations of disorder to be observed experimentally or numerically. In the context of 
Anderson localisation the disorder-averaged moments of wavefunction intensities are readily 
available via standard techniques in the non-linear cr— model framework, see [21 [3] and 
references therein. At the same time extracting typical values of the multifractality exponents 
in that case is a much more challenging task which has been successfully accomplished only 
very recently [3]. In the present paper we would like to concentrate on a different type of 
models where, in contrast, calculating disorder-averaged moments in the full parameter range 
is more difficult, whereas the typical values of IPR exponents are readily accessible. 

Arguably the simplest model with disorder-induced multifractality which attracted 
considerable interest in recent years is the case of a single classical particle subject to a random 
Gaussian potential V{^) logarithmically correlated in space: 



(l/(xi) F(x2)) = -5'ln 



(xi - X2)^ + 
L2 



a <C L, x G (3) 



where we assume |x| < L, and the parameter a stands for a small-scale cutoff. From the point of 
view of equilibrium statistical mechanics the model is characterized by the Gibbs-Boltzmann 
measure P/3(x) = -^^y exp — /3y(x) as a function of the inverse temperature [3 = l/T, and 
the sample size L. The normalization Jj^|^^p^(x)(ix = 1 implies the value of the partition 
function to be given by 



Z{P) = / exp-pV{x)dx. (4) 

J|x|<L 

According to the general discussion, the multifractal structure of the Gibbs-Boltzmann 
measure can be extracted from the knowledge of moments 

Identifying Af ~ (L/a)^ , the Eqs.([5]) and dH) imply the following expression for the typical 
exponents in terms of the appropriately normalized free energy of the system 

Tq = \q\m\m-Qmp\ ^^f^) = - ^^ J^^t^M ■ 

Although the model is well-defined in any A^— dimensional space [Si [7] it is two- 
dimensional situation which attracts most attention, with Eq.(l3|) describing the correlations 
of the Gaussian free field. In particular, for N = 2 the corresponding statistics of Gibbs- 
Boltzmann weights is known to be deeply related to a variety of interesting physical problems, 
ranging from quantum mechanics of Dirac particles in a random magnetic field [HI [9] to 
Liouville model of quantum gravity and theory of self-gravitating particles [IT], see a 
detailed discussion in [6j. Actually, the dimensionality of space plays in many respects only 
secondary role and many (although not all) essential features of the model are expected to 
be A^— independent. The latter point of view is amply supported by the renormalization 
group arguments [6j and by explicit computations in A = oo ^ and A" = 1 [12]. It 
is also worth mentioning that various one-dimensional versions of the problem attracted 
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considerable interest recently in the context of multifractal random walks [TB} I14j . extreme 
value statistics[121 [6l [T5] and quantum gravity-related flQ\ probabilistic questions, see [17] and 
references therein. 



2. Random Energy Model as a toy model for disorder-induced 
multifractality 

A particular extreme "toy model" case of the problem is represented by the famous Random 
Energy Model (REM) by Derrida \18\ I19j. which amounts to replacing random potential by a 
collection of M independent Gaussian variables, after the natural identification M ~ (L/a)^ 
and with the variances scaled with M in the same way as in the logarithmic case: < >= 
2(7^ InM. The only control parameter for the model is 7 = f3'^g'^, and REM is simple enough 
to allow explicit calculation of the free energy |18l I19j . The typical multifractality exponents 
turned out to be given by [9j 



'(?>i 



(7) 



iq-l){l-jq), 0<7< 



1 



q{l-V^r, ^<7<1 • (7) 
0, 7> 1 



The phenomenon of vanishing of the exponents rq>i in the low-temperature phase 7 > 1 
is called freezing and is qualitatively interpreted in terms of the Boltzmann measure being 
essentially localised on a few sites for low enough temperature or strong enough disorder. The 
typical multifractality spectrum corresponding to the above exponents is [9j 

f( ^ I for 7<1 ... 

^^"^ = 1 -^V-4V7«] for 7>1 ' 

where the expression in the first line assumes the range of exponents «_ = (1 — ^77)^ < a < 
1 -I-7, whereas in second line < a < 2y/j. Thinking in terms of the multifractality spectrum 
it is easy to see that the freezing phenomenon at 7 > 1 is related to «_ =0, when the leftmost 
end of the curve f{a) hits the vertical axis precisely at zero level: /(O) = 0. Similarly, the 
change of behaviour of the typical exponent for 7 > 1 /q'^ is induced by dominance of the 
point a_ in the integration over a, Eq.([l]). 

Although obtained in the framework of REM approximation, the above features of the 
typical spectrum are expected to be shared by all the logarithmic model for any > 1 
[9l [6|, which is indeed confirmed by explicit calculations on A^ = 00 [7] and A^ = 1 case 
[12) . In contrast to the case of typical exponents Tq, extracting the "annealed" exponents 
fq from disorder-averaged moments in the logarithmic models poses a serious technical 
challenge. The only systematic attempt in this direction was undertaken for A^ = 2 in the 
framework of mapping to the Liouville model of quantum gravity [10] where it was found that 
fg>i = {q — 1)(1 — 79) for < 7 < 7q = l/(2q — 1) < 1. However, for 7 > 7g the Liouville 
theory was observed to develop unsurmountable singularities and yielded no reliable value of 
the exponents fq. This state of matter clearly calls for reconsidering the problem within the 
general framework of logarithmic models. 

The natural starting point is again the standard REM representing in many respects a 
zero-dimensional limit of the logarithmic models. For such a "toy model" case the disorder- 
averaged moments (IPR's) < Pq > can be evaluated by a well-controlled methods [191 I20j . 
Surprisingly, explicit expressions for the IPR's seem to be available in the literature only in 
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the low-temperature phase 7 > 1 [201 EU [22] . Extending the analysis of [121 1 
M > 1 and q > I: 



we find for 



(P,) = {Z{(3q)/ [Zm' 



M- (9-1) (1-79), 0<7<l/(2g-l) 



M 



¥7 



1+7 
27 



1+7 
27 



r(g) 



r 



r(g)r 1 



1 

V7 



2? 



7 > 1 



^ <7< 1 



(9) 



where we included in the last line the well-known low-temperature results of |20l 1211 [22] . with 
r(x) standing for the Euler gamma- function. The fact that for 7 > 1 the moments remain finite 
in the limit M — > 00 reflects the quasi-localised nature of the Boltzmann-Gibbs measure in 
the low-temperature phase. As a consequence, the "annealed" multifractal exponents remain 
frozen: fq = 0. In the high-temperature phase the exponents fq are non-vanishing and in the 
range < 7 < ^ typical and annealed exponents coincide. The annealed exponents actually 
keep that common value up to 7 = 7^ = l/(2g — 1). Both the value fq for 7 < 7g and the value 
of the threshold 7g are in full agreement with the Liouville model analysis [lOj. Finally, in the 
range 7^ < 7 < 1 the annealed exponents change drastically by acquiring the g— independent 



(1-7)^ 



The behaviour of typical and 



value which relates to 7 in a non-polynomial way fq — — 
annealed multifractal exponents for various values of the control parameter 7 is summarized 
in the diagram, and is further discussed below. 



Tn 



(g-l)(l-7g) 

simple 
scaling 



9(1 - ViY 

- dominance 





{q-m-iq) 
simple scaling 



l/(2g 

— r 





freezing 



7 



1 - 7)727 
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freezins 



Figure 1. Regimes of behaviour for two types of multifractal exponents for the 
Random Energy Model: typical exponent Tq vs. annealed fq at different values of 
effective disorder 7. 



Using Eq.([9]) one can restore the corresponding mean density of multifractal exponents: 

/ M \ 



< p{a) >= ( ^ 5(a - ai) ) w C(M, a) M-^(") 
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where the "annealed" multifractality spectrum turned out to be given by 

, f l-^[a-(l+7)]', for 0<7<1 and < a < 1 + 7, , . 
\ -^[a^-4^a] for 7>1 and < a < 2^ ' 

These expressions show that the disorder- averaged (or "annealed") spectrum is precisely the 
same as the typical spectrum corresponding to Eq.Q, with the only essential difference 
being that annealed spectrum in the first line of Eq. ijlOp becomes negative in the range 
< a < a_ = (1 — ■^/7)^- These values of a correspond to exponentially rare events in 
full agreement with our earlier discussion and the picture developed in [21 [5]. The value 
a = is the lowest possible value of the exponent a, and for this reason is frequently 
called the "termination point" of the (disorder-averaged) multifractality spectrump]. Further 
substituting the corresponding < p[a) > to the integral over cfl, see Eq.©, we find that in 
the range ^q-^i < 7 < 1 the integral for (Pqyi) is actually dominated in the limit M 1 by 
the vicinity of the lower integration limit a = rather than by the stationary point of the 
integrand. It is in this way that f, acquires the q-independent value cf. Eq.®, as was 

indeed anticipated in [2j on heuristic grounds. We see that the actual mechanism behind the 
observed drastic change of the annealed exponents is very close to one forcing the exponent 
freezing at 7 > 1, and in a sense is the precursor of the latter behaviour. For this reason it 
is natural to suggest to call this phenomenon pre-freezing. We shall also see later on that the 
replica approach gives additional support to associating the observed behaviour with a partial 
freezing of certain kind. 



3. Annealed multifractality exponents for the circular logarithmic model 

After the detailed understanding of the toy REM limit it is natural to try to extract annealed 
exponents for a more realistic case of one-dimensional model with logarithmic correlations. 
The most promising instance of ID system of that type is arguably the "circular logarithmic 
model" (CLM) introduced in [12j . To define it, consider the lattice of M points positioned 
equidistantly at the circumference of a circle of a radius a <^ R < L in the standard two- 
dimensional free field with correlations Eq.([3|). The angular coordinates of the points are 
given by 6*^ = ^k, k = 1,2,... ,M. Then it is easy to to see that the values of the free field 
Vi associated with those points are characterized by covariances 

Cki = {VkVi) = -g' In |4sin2 + 2g^ In {L/R), k / I (11) 

The first term in (jlip defines precisely CLM as described in p2], with constant second term 
playing a somewhat trivial role and thus omitted henceforth. It turns out that the consistency 
of the procedure requires to choose the variance (V^) = V"^ to satisfy V"^ > 2g'^liiM, and 
following [12] we choose simply V'^ = 2g'^ InM. The partition function for our model is defined 
in the standard way through Zp = e~^'^' , with the goal to evaluate the IPR moments 

Pq = {Z{j3q)/ [Z{(3)]'^) in the limit M S> 1 for various values of the control parameter 7 = (3'^g'^. 

In the rest of the paper we restrict our explicit calculations by the simplest representative 
case q = 2, and employ a variant of the replica trick: 

= {zm/ [Zm') = lim (z(2/3) [ZiP)r~') , (12) 

II It is worth mentioning that for REM one can extract not only the multifractality spectrum, but 
also the expressions for the pre-exponential factors C{M, a) which goes beyond the standard precision 
of the multifractality analysis. For example, in the high-temperature phase < 7 < 1 one finds 
C{M,a) - iyCT (i_ M--)-^. 
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implying a kind of continuation from integer values n > 2 to n = 0. The key point is that 
for integer n > 2 the disorder averaging in (jl2p can be performed by the method described in 
detail in [12], which gives, in particular 

(z(2/3)[Z(/3)]'=)|Af»i = Mi+'=+^('=+^)jfc(7), 0<A:<i-3 (13) 

where Jk{l) is the Dyson-Morris Integral [23] given by a product of gamma- functions: 

Ml) = 7:^ de,... d9,l{ le'''^ - e''^ e^'' (14) 

(2vr)'=yo Jo fi' 

1 r[l-7(fc + 2)]r[l-7(A: + 3)] 

[r(i-7)]fc-ir[i-7(A: + i)]r(i-27)r(i-37) ■ ^ ' 

Performing in the above expression the "naive" replica limit k ^ —2 we arrive at the following 
expression for the IPR: 

^ r(i + 7)r(i - 27)r(i - 37) ^ ' 

We see that the value of the annealed multifractality exponent f2 = 1 — 27 which emerges from 
our calculation coincides with the "simple scaling" value Tg = {q — 1)(1 — qj), q = 2 discussed 
by us earlier. Moreover, the consistency of the above procedure obviously requires < 7 < 1/3, 
with the upper limit being precisely the threshold 7^=2 = 1/3 of validity of the above "simple 
scaling" regime. We conclude that a simple-minded replica limit could be employed to produce 
meaningful results only as long as the pre-freezing phenomenon responsible for the change of 
multifractality exponent is not operative. 



4. Infinite-dimensional limit: prefreezing via replica symmetry breaking 

To get some progress in understanding of the mechanisms behind the failure of the simple 
scaling in replica approach, we turn from now on to another exactly solvable limit of the 
logarithmic model, that is to the infinite-dimensional case. The free energy, hence the typical 
multifractality spectrum, was calculated in [7j in the framework of the replica trick, and 
very recently confirmed by rigorous mathematical methods [24j. The system was found to 
display the REM- type freezing transition at 7 = P'^g'^ = 1, with the low-temperature phase 
7 > 1 described by the standard one-step replica symmetry breaking pattern. The meaning 
of the freezing could be elucidated by invoking the probability for two independent particles 
distributed in such random potential according to the Boltzmann-Gibbs measure to end up at 
a distance of order of the small cutoff scale a"^. The probability was found [7] to tend to zero 
in the thermodynamic limit L/a — > 00 everywhere in the high-temperature phase < 7 < 1, 
confirming the particle delocalization over the sample. In contrast, in the lower temperature 
phase 7 > 1 two particles have a finite probability to be trapped at the small-scale distance 
even in the infinite sample. 

Our starting point here is again the identity (jl2p . Employing it one can easily perform 
the disorder average for any integer number of replica n > 2. After appropriate rescaling of the 
coupling constant g — > g^/N and length scales L — > lV^N and a — > a^/2N the manipulations 
similar to those described in detail in [7j yield a convenient representation for the IPR in terms 
of an integral over a positive definite matrix Q of the size (n — 1) x (n — 1) with entries q^^u- 
We have 

(P2) = limCjv,n(a)i^"' / (detQ)~"/2 g-Ar*„{Q) (^7) 
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where 

^n(Q) = -^ln(detQ)+7 In 



2 

l<fM<u<n-2 



+ 27 ^In 



1 , , 2' 

2\1fi,fJi + 1n-l,n-l) — Qti,n-1 + 0, 



7 = /5V, (18) 



and the integration domain Dq in the above expression is simply Dq = {Q > 0, q^^^ < 
L^, ij, = 1, ... n — 1}. The proportionality constant CN,n{cL) is also explicitly known, but its 
value is inessential for the subsequent calculation. 

The shape of the integrand in p!7|) is suggestive of application of the saddle-point method 
for evaluation of the integral in the large- limit. The corresponding saddle-point equations 
for the entries of the matrix Q amount to d^/dq^^i, = for any choice of the indices 
l<;U<i^<(n — 1). A closer inspection of the replica limit n — > reveals however that 
solutions to the saddle-point equations do not actually exist unless one fixes all the diagonal 
entries g^.^, 1 < /i < (n — 1) of the matrix Q to the boundary of the integration domain by 
setting q^^^ = L^, and excluding them from the variational procedure ( cf. a similar result in 
[7]). The remaining off-diagonal entries should be found from the system of equations: 

[Q''l.- + ^ L^-,'. + a' -°- l£M<='<"-2 (19) 

[Q"l.,.-i + ^^ L^-J..+„. =°. ISA'S "-2 (20) 

The apparently non-equivalent roles played by the special replica index = n — 1 and 
the rest of n — 2 replicas in the above equations are inherited from the structure of IPR 
representation in the replica method, Eq. (|12|) . The only ansatz for Q respecting the full 
permutation symmetry between the (n — 2) equivalent replicas is given by 

Qr... = ( ) , Q^°\. . = {L^ - qo)S^,u + go, 1 < /i, 1/ < n - 2 (21) 



where Q^^^^ > is of the size (n — 2) x (n — 2), and the (n — 2)— component vector v is of 
the form v = 11(1,. ..,1). Two parameters qo and v characterising such a replica-symmetric 
solution satisfy in the limit n ^ the system of equations 

qo - v'^/L'^ 7 



(L2 - go)(L2 - 3<?o + 2v^/L^) (L2 - q^ + a^) 

V 27 
L2(L2-3go + 2t;7L2) " (L2-^ + a2) ' 



(22) 
(23) 



The only solution of such a system existing in the high-temperature phase in the 
thermodynamic limit L ^ a has the following form 

^0 = ^J^^L' + O{a% v = -^L' + 0{a'). (24) 
(I + 7) 1 + 7 

The condition qo < which is readily seen to be satisfied everywhere in the high-temperature 
phase < 7 < 1 ensures that Q^''^ > 0. 

Substituting these values back to the relation (fTH]) and using the identity 

det Q,.,. = L\L^ - go)"-' [L^ - v^/L" + {n - 2,){qo - v^/L^)] (25) 
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we find to the leading order ^n^oiQr.s.) = (1 ~ ^7) InL whicli in turn yields the asymptotic 
behaviour for the averaged inverse participation ratio (P2) oc e-^*"-*o oc M-(i-2t) where 
we have identified M ~ . We arrive to an important conclusion of the replica symmetric 
ansatz Eq. (|2ip implying necessarily the "simple scaling" result f2 = 1 — 27 for the annealed 
multifractality exponent, see the discussion after (jl2p . As is clear, the latter value cannot have 
any meaning at least for 7 > 1/2, which in turn implies that the replica-symmetric solution 
of the problem cannot be valid in the whole high-temperature phase < 7 < 1. The only way 
out is therefore to look for an appropriate scheme of replica symmetry breaking which occurs 
already in the high-temperature phase. 

To get a guiding idea in our search for a solution to the saddle-point equations (jl9l20p 
which goes beyond the replica symmetric ansatz, eq. (j2ip . it is useful to recall that the singled- 
out replica indexed with ly = n — 1 has its origin in representing the factor Z{2[3) in the 
averaged inverse participation ratio, cf. Eq. (fT2|) . As such that special replica is effectively 
"colder" then all other equivalent replica indexed with 1,2, ... ,n — 2 which originated from 
the factors Z{0). It is then natural to suppose that by interacting with other replica the colder 
one could eventually "pre- freeze" a certain group of replica around it. The manifestations of 
the freezing mechanism within the replica approach was discussed for the present model in 
detail in [7]. Employing it one should expect such a "pre-freezing" to manifest itself via 
existence of a diagonal block inside the matrix Q having in the thermodynamic limit L ^ a 
all entries equal , up to the leading order, to L^. 

The simplest ansatz for the matrix Q compatible with a possibility of such pre-frozen 
group of m replicas, < m < n — 2, would be of the following structure: 

O -(Tv 0^2)^ S\ Q^'\,u = iL'-qi)S,,u + qi, 1 < f^, 1^ < n - m - 2 
Hr.s.b — ^S-t- ' V2 , ^(2) /r2 \x , i ^ ^ i^"-' 

\ ^2 j V,^ = ~ Q2)S^,,u + 92, l<n,i'<m 

where all entries of the matrix E of the size m x [n — m — 2) are equal to 1, the 
(n — m — 2)— component vector vi has all components equal to vi and m— component vector 
V2 is of the same structure: V2 = ^2(1, •••)!)• The values of five parameters qi,q2, q3,vi,V2 
for a given size m of the pre-frozen block are to be found from the saddle-point conditions 
(jl9p . (j20p in the replica limit n ^ 0. Finally, the parameter m satisfying in the replica limit 
inequality —2 < m < should itself be chosen as to extremize the resulting ^(Q), eq. (|18p . 
which is evaluated with the help of the identity 

det Qr.s.b. = L\L^ - gi)"-'"-^ [L^ -vf/L^ + {n-m- 3){qi - vf/L^)] 
X - q^r-^ - vj/L^ + {m- l)iq2 - vj/L^) -rm{qs - viV2/L^)] (27) 

n — m — 2 

where r = -— n—^ ; o-tt^ ■ 

L2 -vl/L^ + {n-m - 3){qi - vf/L^) 

Substituting the Ansatz (p6]) into the saddle-point equations (fT9]) . (pOj) one arrives in the 
limit n = after straightforward but lengthy algebraic manipulations to the set of equations: 

7 ' r +7"^ 72 I 2 + 72 , 2 (^8) 

" — 53 + — f 1 + 

+ irn 2 +^7 ,2 2 (29) 

^ Q2 + o L'' — V2 + a 

7(m + 2) , +27 , (30) 



Qi 


L2- 










q2 


Q3 





L2_ 


qi{m + 3) 


L2 


- gi + a2 


L2_ 


qi{m + 3) 


L2 


- qs + a"^ 


L2 + 


■ q2{m - 1) 



L2 — qi — q^ + a2 — qi + — vi + 
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-i—r9. — — - 7 ("i + 2) — ——0+^1 TO , o (31) 



— q2 L"^ — q2 + L"^ — + L"^ — V2 + 

^^1 = ^7^^ r^+^7"i72 (32) 

— f 1 + — i;2 + a 

l = ZjL^ + 27(m + 2) a '^Tm 2 ' (^3) 

-L^ — f 1 + — i'2 + 

where we found it convenient to use the parameter D = [Q^"'^]„_i^„_i. 

It is immediate to check that the above system of equations always has the rephca 
symmetric type of solution such that = 92 = 93 = qo and vi = V2 = v. Indeed, such 
a substitution results in m dropping out from the equations, and go and v satisfying equations 
(j22p . At the same time there exists another type of solution which explicitly depends on m 
and is given in the limit L ^ a by 

qi = ^^l^Y [1 + 7("i' + 3m + 3)] + ©(a^), q2 = + 0{a^) 

q^ = --^{m + 2)L^ + 0{a^), vi = --^{m + 2)L^ + 0{a^), V2 = L^ + 0{a^). (34) 
i + 7 i + 7 

We indeed see that all the off-diagonal entries inside the block i are such that 



V2 L\ 

the (squared) "distances" d = q2 — and d = v — are of the order of the cut-off scale 
a?. The latter feature is precisely the property of being frozen, see [7j for a more detailed 
discussion. It is also worth noting that though in the limit m — > in (j34p qi and vi tend to 
the replica symmetric values, they remain different from 92, 93 and V2, respectively, so that 
even in the limit m — > the new solution does not tend to the replica-symmetric one. 

Substituting now the solution (j34p into (j27p one finds after straightforward algebraic 
manipulations | Indet Qr.s.fe. = — (1 + m)lnL + O(lna), which after more algebra results in 
the value of the functional <I>„^o(Q) foi" a given value of the parameter m: 

^n^o(Qr.s.6) = In L (l + m - -i[n? + 3m + 2]) + 0(ln a) (35) 

Extremizing this functional with respect to m G [—2,0] one immediately finds that for 
1/3 < 7 < 1 the extremum is at an internal point of the interval m = = whereas 
for < 7 < 1/3 the extremum is at the boundary point m = 0. Substituting those values to 
the functional leads in the thermodynamic limit to value <I>„^o(Q) given in such a scheme to 
the leading order by 

f 1 - 27, < 7 < 1/3 

\ 1/3<7<1 

The above result immediately implies the values for the corresponding multifractality exponent 
7^2 coinciding everywhere in the high-temperature phase < 7 < 1 with that exemplified by 
the standard REM, cf. dH). 

To complete the picture, one can perform the standard de-Almeida-Thouless-like 
analysis [25] and verify that the replica-symmetric solution Eq. (j24p does not show any local 
instability at the point of the pre-freezing transition 7^1/3 — (although it does become 
unstable for a certain value 7 = 'jinst £ (1/3, 1/2))- Absence of local instability is consistent 
with the mentioned observation that in spite of m — > when approaching the pre-freezing point 
the solution with broken replica-symmetry ()34p remains formally different from the replica- 
symmetric solution (|24p . Above the pre-freezing temperature, that is for 7 < 1/3, the two 
solutions share the same value of the functional ^\extr, and therefore produce exactly the same 
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simple scaling for the multifractality exponent f2 . In the pre-freezing domain 1/3 < 7 < 1 the 
solution with broken replica symmetry wins. 

In conclusion, we have performed systematic analysis of the multifractality exponents 
extracted from the averaged moments of the Boltzmann-Gibbs measure generated by 
logarithmically correlated random potentials. In particular, using zero-dimensional and 
infinite-dimesnional versions of the model we have identified a pattern of the replica symmetry 
breaking responsible for the abrupt change ("pre-freezing") of those exponents in the high- 
temperature phase. Implementing such a pattern in explicit calculations for one- and two- 
dimensional versions of the models remains a challenging open problem. 

Acknowledgements. This work was supported by Leverhulme Research Fellowship project "A 
single particle in random energy landscapes" . Major part of the research has been completed 
during the author's participation in the programme on Mathematics and Physics of Anderson 
Localisation at the Newton Institute(Cambridge) whose support and hospitality is gratefully 
acknowledged. The author is grateful to A. Mirlin for stimulating initial discussions on the 
problem of annealed multifractality exponents, and to C. Mudry, P. Le Doussal and M. 
Skvortsov for encouraging interaction at various stages of the work. 



[1] Paladin G and Vulpiani A Phys. Rep. 156 (1987) 147 

[2] Evers F and Mirlin AD, Rev. Mod. Phys. 80 (2008) 1355 

[3] Foster MS, Ryu S, and Ludwig AWW (2009), e-preprint arXiv:0901. 02841 

[4] Monthus C and Garel T Phys. Rev. E 75 (2007), Art. No. 051122 

[5] Mandelbrot B Physica A 163 (1990), 306; Chabra AB and Sreenivasan KR Phys. Rev. A 43 

(1990) 1114 ; Halsey TC, Honda K, and Duplantier B J. Stat. Phys. 85 (1996) 681 
[6] Carpentier D and Le Doussal P Phys. Rev. E 63 (2001), 026110 

[7] Fyodorov YV and Bouchaud JP J. Phys. A: Math.Theor 41 (2008) 324009 (25pp); Fyodorov YV 

and Sommers H-J Nucl. Phys. B [FS] 764 (2007), 128 
[8] Ludwig A, Fisher MPA, Shankar R, and Grinstein G Phys. Rev. B 50 (1994), 7526 
[9] Chamon C, Mudry C and Wen X-G Phys. Rev. Lett. 77 (1996) 4194; Castillo H E, Chamon C C, 
Fradkin E, Goldbart P M, and Mudry C Phys. Rev. B 56 (1997) 10668 
[10] Kogan II, Mudry C, and Tsvelik AM Phys. Rev. Lett. 77 (1996) 707 
[11] Abdalla E and Tabar M R R Phys. Lett. B 440 (1998) 339 

[12] Fyodorov YV and Bouchaud JP J. Phys. A: Math.Theor 41 (2008) 372001 (12pp) 

[13] Muzy J-F, Delour J, and Bacry E Eur. Phys. J. B 17 (2000) 537 Bacry E, Delour J, and Muzy 

J-F Phys. Rev. E 64 (2001) 026103 
[14] Ostrovsky D Lett. Math. Phys 83 (2008) 265 
[15] Fyodorov YV, Le Doussal P, and Rosso A, under preparation. 
[16] Duplantier B and Sheffield S (2008), e-preprint arXiv:0808.1560 [math.PR] 
[17] Rhodes R and Vargas V (2008) e-preprint arXiv:0807. 1036 math.PR] 
[18] Derrida B Phys. Rev. B 24 (1981) 2613 
[19] Gardner E and Derrida B J. Phys. A 22 (1989) 1975 
[20] Derrida B and Tolouse T J Physique Lett 46 (1985) L223 
[21] Derrida B Physica D 107 (1997) 186 

[22] Bouchaud J-P and Mezard M J. Phys. A: Math. Gen. 30 (1997) 7997 

[23] Forrester P J and Warnaar S O , Bull. Amer. Math. Soc. (N.S.) 45 (2008), no. 4, 489-534 

[24] Klimovsky A (2009), "Parisi landscapes in high-dimensional Euclidean spaces" , talk at the workhop 

"Mathematical Models from Physics and Biology", April 2009, Bonn, Germany 
[25] de Almeida JRL and Thouless DJ J. Phys. A 11 (1978), 983 



